Large Fluctuations, Classical Activation, Quantum Tunneling, and Phase Transitions 
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We study two broad classes of physically dissimilar problems, each corresponding to stochastically 
driven escape from a potential well. The first class, often used to model noise-induced order param- 
eter reversal, comprises Ginzburg-Landau-type field theories defined on finite intervals, perturbed 
by thermal or other classical spatiotemporal noise. The second class comprises systems in which a 
single degree of freedom is perturbed by both thermal and quantum noise. Each class possesses a 
transition in its escape behavior, at a critical value of interval length and temperature, respectively. 
It is shown that there exists a mapping from one class of problems to the other, and that their 
respective transitions can be understood within a unified theoretical context. We consider two ap- 
plications within the first class: thermally induced breakup of monovalent metallic nanowires, and 
stochastic reversal of magnetization in thin ferromagnetic annuli. Finally, we explore the depth of 
the analogy between the two classes of problems, and discuss to what extent each case exhibits the 
characteristic signs of critical behavior at a sharp second-order phase transition. 

PACS numbers: 05.40.-a, 02.50.Ey, 03.50.Z, 68.65.La, 75.60.Jk 
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I. INTRODUCTION 



Noise-induced escape from a locally stable state gov- 
erns a wide range of physical phenomena yj . In spatially 
extended classical systems, where the noise is typically, 
but not necessarily, of thermal origin, these phenomena 
include homogeneous nucleation of one phase inside an- 
other Q , micromagnetic domain reversal 0, 0] , pattern 
formation in non-equilibrium systems Q , and many oth- 
ers. The concepts and techniques employed are formally 
identical to those used in tunneling problems in quantum 
field theories, including among others the 'decay of the 
false vacuum' [|| , anomalous particle production 01 , and 
macroscopic quantum tunneling in a dc SQUID 



Recently a new kind of 'phase transition', in the clas- 
sical activation behavior of spatially extended systems 
perturbed by weak spatiotemporal noise, has been uncov- 
ered and analyzed |3, |lfj, |ll| . The existence of a similar 
crossover, from thermal activation to quantum tunneling 
in systems as simple as a single degree of freedom, has 
long been known [llBElM EKE GUM 
Here we will review the stochastic escape behaviors of 
both types of systems, explore the depth of the anal- 
ogy between them, discuss to what extent the crossover 
can be thought of as a type of second-order phase tran- 
sition, and show how the purely classical transition, oc- 
curring as system length is varied, can shed light on the 
crossover from classical activation to quantum tunneling 
as temperature is lowered. But let's begin with a purely 
classical problem. 



II. MODEL 

Consider an extended system described by a classical 
field 4>(z,t) defined on the spatial interval [-L/2, L/2], 
subject to the potential V((p) and perturbed by spa- 
tiotemporal white noise. Its time evolution is governed 
by the stochastic Ginzburg-Landau equation 



(1) 



dtt = d zz - d$V(<f>) + V2T£(z, t) , 



where all dimensional quantities have been scaled out. 
The first term on the RHS arises from a field 'stiff- 
ness', i.e., an energy penalty for spatial variation of the 
field. The system is stochastically perturbed by addi- 
tive white noise £(z,t), satisfying (£(zi, ti)£(%2, £2)) = 
5(zi — Z2)S(ti — £2)- We consider only weak noise, i.e., 
the noise magnitude T is small compared to all other en- 
ergy scales in the problem (formally, our analysis will be 
asymptotically valid in the T — * limit). The frictional 
coefficient has been set to one, so when T is temperature 
Eq. 0) obeys the fluctuation-dissipation relation. 

The deterministic, or zero-noise, dynamics can be writ- 
ten as the variation of an action H with the field 6: 



<p = sn/S(j), 



(2) 



and because the dynamics is gradient, the action TL is 
simply an energy functional 



n\d>] = 



L/2 



dz 



L/2 



(3) 



We will mostly consider relatively simple potentials, in 
particular the symmetric bistable quartic potential 



1 



1 



14W = -/+/, (4) 
and the asymmetric monostable cubic potential 
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vm = \<P -\f - (5) 
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FIG. 1: Potentials used in discussion: (a) bistable quartic 
(Eq. Q) and (b) monostable cubic (Eq. JSJ). 



Both potentials are shown in Fig. ^ 

In the weak-noise (T — > 0) limit, the classical activa- 
tion rate for a transition out of a stable well is 



r-r exp(-A J B/T). 



(6) 



where AE is the activation barrier (to be defined below) 
and Fo the rate prefactor, determined by fluctuations 
about the most probable escape path. When the top 
of the barrier is locally quadratic, the prefactor r is 
independent of temperature. In such circumstances the 
escape rate © is said to be of the Arrhenius-van't Hoff 
(or often simply Arrhenius) form. 



III. THE INFINITE-DOMAIN CASE 

Let's start by considering the L = oo case, which was 
fully worked out for the classical nucleation problem by 
Langer and for the quantum tunneling problem by 
Callan and Coleman 0. (See also Schulman |23j for a 
good pedagogic treatment.) For illustrative purposes, we 
will use the bistable quartic potential l@J). 

The states of interest — i.e., the stable and saddle 
field configurations — are time-independent solutions of 
the zero-noise Ginzburg-Landau equation <fi = —67i/S<p. 
They are therefore extrema of the action so in the current 
case satisfy the nonlinear differential equation 

r = -cp + 4?. (7) 

Its spatially uniform solutions are just the stable states 
4>{z) — ±1 and the local maximum cf> = 0. But our main 
interest is in the soliton-like pair (or 'bounce' (f|) 



= ± tanh 



(z-z Q )/V2 



(8) 



where zq is a constant whose presence denotes the fact 
that the 'domain wall', i.e. the spatially varying piece 
of JSJl that separates the two stable states, can nucleate 
anywhere on the line. 




FIG. 2: The 'bounce' described by Eq. (JSJ with zo 
clarity, only one of the symmetric pair is shown. 



0. For 



Assuming for the moment that ijBJ is the saddle con- 
figuration (this will be justified below), the activation 
energy AE = -E[0 sa ddic] - instable] in the Kramers rate 
formula © can be computed, giving 



c °° 

AE= I dz 



:(d z sa ddlc) 2 + V((j> sadd l e ) 



dz 



stable ) 



stable J 



2V2 



(9) 



which is essentially the energy of a single domain wall. 

In order to compute the prefactor Tq of the Kramers 
rate 10, we need to examine fluctuations about the op- 
timal escape path, in particular in the vicinity of the 
stable and saddle field configurations. The procedure for 
doing this is described in detail elsewhere (see, for exam- 
ple, [IBl23)l23)> an d w iH be simply summarized here. 

Let ip s denote the stable state (in this case, the uni- 
form ±1 state), and let <pt denote the transition (saddle) 
state (here given by Eq. ijSJ). Consider a small perturba- 
tion r\ about the stable state, i.e., ip = ip s + Then to 
leading order i\ = —A s rj, where A s is the linearized zero- 
noise dynamical operator governing the time evolution 
of fluctuations about tp s . Similarly A t is the linearized 
zero-noise dynamical operator for (p t . 

We next diagonalize the linear time-evolution opera- 
tors, by decomposing fluctuations about the stable and 
transition states into normal modes, which are eigenfunc- 
tions rji of the cooresponding operators: 



(10) 



where b — s, t. An eigenfunction with positive eigenvalue 
Xi > is a stable mode; one with Aj < is unstable. 
A stable (or metastable) state therefore has all > 0; 
a saddle state has a single A; < 0. Its corresponding 
eigenfunction denotes the unstable direction (in function 



3 



space) in the vicinity of the saddle, leading either back to 
the initial stable configuration or out of the well entirely. 
With the potential (J2J, the eigenvalue equation becomes 



fi = -A b r) = - [-d 2 /dz 2 + (-1 + 3<fi)] n ■ 



(11) 



In most cases the barrier is locally quadratic: all eigen- 
values are nonzero and we're left with an infinite set of 
decoupled quadratic fluctuations about the stable and 
saddle states. Then 



To = 



2tt' 



det A 6 



det A t 



|A*, | 



(12) 



where At ; o is the only negative eigenvalue of At. In gen- 
eral, the determinants in i|12[) separately diverge: they 
are products of an infinite number of eigenvalues with 
magnitude greater than one. However, their ratio, which 
can be interpreted as the limit of a product of individual 
eigenvalue quotients, is finite. 

There is a technical difficulty that has not yet been 
addressed, namely the existence of a soft collective mode 
corresponding to the arbitrariness of zq in JHJ: the do- 
main wall (or 'instanton') can nucleate anywhere. The re- 
sulting translational symmetry implies a zero-eigenvalue 
mode. Its removal can be achieved with the McKane- 
Tarlie regularization procedure [2(j (see also [22j]) for 
functional determinants. This will not be discussed fur- 
ther here, except to note that its presence results, among 
other things, in a non-Arrhenius (noise-dependent) pref- 
actor that scales with the length. The overall result for 
the prefactor per unit length is then 



r„M^)( 



(13) 



where the first term on the RHS follows from the com- 
putation of Eq. (|12[1 with the zero eigenvalue removed 
via the McKane-Tarlie procedure, and the second term 
(divided through by length L) gives the contribution of 
the zero eigenvalue, i.e., the effect of the translational 
symmetry of the bounce. 



IV. THE FINITE-DOMAIN CASE 

The differences between the infinite and finite-line 
cases are not only quantitative; there are important 
qualitative differences that are also relevant, in an en- 
tirely different context, to the classical *-y quantum 
crossover (Sec. IVIf) . The most striking of these differ- 
ences is a sharp change in activation behavior as interval 
length is varied; moreover, this change exhibits charac- 
teristics of a second-order phase transition jjj, [l(J, [flj , but 
only in a strictly asymptotic sense. 

Because we are now working on an interval of finite 
length, we need to specify boundary conditions. With 
the exception of the zero-energy mode that accompanies 
only translation-invariant (such as periodic) boundary 



conditions, different choices of boundary conditions lead 
only to minor quantitative differences. To avoid the (mi- 
nor) complication of the zero mode altogether, we will 
employ Neumann boundary conditions: d4>/dL\_ L / 2 — 
dcf)/dL\ L /2 — 0. With this choice Stable = ±1, as before. 
We will continue to use the notation of Sec. lIIII where 4> s 
(</>t) refers to the stable (transition) state. 

For the symmetric <p potential with Neumann bound- 
ary conditions, the change in activation behavior arises 
from a bifurcation of the transition state, from a uniform 
configuration below a critical length L c to a pair of de- 
generate, spatially varying 'bounce' configurations above 
L c . More precisely, the saddle states are 







when L < L c and 



± 



2m 



1 



■sn( 



to) , 



(14) 



(15) 



when L > L c , where sn(- | to) is the Jacobi elliptic 
sn function with parameter < m < 1. Its quarter- 
period is given by K(to), the complete elliptic integral 
of the first kind |28j| . which is a monotonically increas- 
ing function of to. As m — * + , K(to) decreases to 7r/2, 
and sn(- | to) — * sin(-). In this limit the saddle state 
smoothly degenerates to the (f> — configuration. As 
to — > I — , the quarter-period increases to infinity (with 
a logarithmic divergence), and sn(- | m) — » tanh(-), 
the (nonperiodic) single-kink sigmoidal function. The 
Langer/Callan-Coleman bounce solution (jHJ) is thereby 
recovered as L — > oo. 

The value of to in l|15[) is determined by the interval 
length L and the Neumann boundary conditions, which 
require that 



L/Vm+1 = 2K( to ) . 



(16) 



The critical length is determined by (|f 6[) when to = 0; 
that is, L c — 7r. As previously noted, to — > 1 corre- 
sponds to L — ► oo, and the activation energy smoothly 
approaches the asymptotic value of 2y2/3. The transi- 
tion state for an intermediate value of to, corresponding 
to L = 10, is shown in Fig. [3] 

The activation energy AE can be solved in closed form 
for all L > L c (below L C1 it is simply L/4): 



1 



AE = —, 7^7 4(1+to)E(to)--(1-to)(3to+5)K(to) 

3(1 + to) j ' a L 2 

(17) 

with EJro) the complete elliptic integral of the second 
kind [28j. The activation energy as a function of L is 
shown in Fig. Note that the curve of AE vs. L and 
its first derivative are both continuous at L c ; the sec- 
ond derivative, however, is discontinuous, as might be 
expected of a second-order-like phase transition. 

A more profound manifestation of critical behavior at 
L c is exhibited by the rate prefactor To. When L < L c = 
7r, both the stable and saddle states are spatially uniform: 



4 




FIG. 3: The transition state (j>t(z) for L — 10 (corresponding 
to m = 0.986) described by Eq. I|15|l . As in Fig. EI on ly one 
of the symmetric pair is shown. 



AE 




10 



12 



14 



FIG. 4: The activation energy AE as a function of the in- 
terval length L, for the potential given by Eq. Q (with all 
coefficients set equal to one) and Neumann boundary condi- 
tions. The dashed line indicates the critical interval length 
L c — 7T at which the saddle state bifurcation takes place. 



^stable = ±1 and sac idic = 0. This greatly simplifies the 
computation of the associated eigenvalues. 

It is immaterial which of the two (degenerate) stable 
states is used; the eigenvalue spectrum is the same at 
both because of the symmetry under (j) i— > — <fi. Lineariz- 
ing around either stable state yields the operator 



-d 2 /dz 2 



(18) 



and similarly 

A t = -d 2 /dz 2 - 1 . (19) 
The eigenvalue spectrum of A s with Neumann boundary 



conditions is 

•^s,n 2 



L 2 



n = 0,1,2. 



The eigenvalue spectrum of A t is similarly 

2 2 

\ t =-l + —— n = 0,1,2.. 



(20) 



(21) 



As required, all eigenvalues of A s are positive, while A t 
has a single negative eigenvalue A^o = —1- Its eigen- 
function, which is spatially uniform, is the direction in 
configuration space along which the optimal escape tra- 
jectory approaches 4>t- 

Putting everything together, we find the Neumann- 
case rate prefactor when L < L c to be 



Z7T 



n 



n=0\ 



IT 



n 



n=0 



IT 



h(V2L) 



' sin 



2 3 / 4 7 



sinL 



(22) 



As L -> L~{= 7I-), T - (L c - L)- 1 ' 2 . This divergence 
has a simple physical interpretation: the optimal escape 
trajectory becomes transversally unstable, in the direc- 
tion defined by the eigenmode rji, as the critical length 
is approached. Mathematically the divergence is caused 
by A t: i — > + as L — > L~ . 

When L > L c , there are two transition states, namely 
the nonuniform bounce configurations ±0 t given by 11511 . 
The associated linearized evolution operator, computed 
from (|llfl . is 



' dz 2 



1 



6m 



1 



\Jm, + 1 



(23) 



Calculation of the associated determinant quotient, 
and the sing le unstable eigenvalue Xt,o, is described in 
detail in [l(|, to which the interested reader is referred. 
The Neumann-case rate prefactor when L > L c is 



r.-I 



1 - 



-\J m 2 



1 



sm 



h(V2L) 



72 1(1 - m)K(m) - (l + m)E(m)| 



(24) 



As m -> 0+ (L L+), T - {L-Lc)' 1 ' 2 . The prefactor 
over the entire range of L is shown in Fig. |3J 

The divergence at L c is striking, but requires interpre- 
tation. We defer further discussion to Sec. lVIII and turn 
now to physical applications of the methods and results 
presented in this section. 



V. TWO APPLICATIONS 

We have established a transition in activation behavior 
as system length is varied; but is it, and more generally 
the theory presented above, relevant to actual systems? 
In this section we briefly discuss two applications. 
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FIG. 5: Rate prefactor To for the quartic potential with Neu- 
mann boundary conditions, showing the power-law divergence 
of the prefactor as L — > L^r. 



A. Lifetimes of monovalent metallic nanowires 

Metallic nanowires are cylindrically-shaped incom- 
pressible electron fluids with diameters of order tens of 
atoms and with lengths hundreds to thousands of atoms. 
They are stabilized by quantum shell effects [23 . 18(1 l3rj 
but at nonzero temperatures are only metastable, with 
breakup probably due to thermal fluctuations |32 . l33 . l3^ |. 
We have proposed [35| a self-consistent continuum ap- 
proach to studying the lifetimes of monovalent metallic 
nanowires, with a large-deviation-induced 'collapse' mod- 
elled through a stochastic Ginzburg-Landau field theory, 
of the kind discussed above. Our theory provides good 
quantitative agreement with available data on nanowire 
lifetimes, and accounts for the observed difference in sta- 
bility between alkali and noble metal nanowires. 

We treat a nanowire as a cylinder of length L and ra- 
dius R(z) = R + 4>{z), with z £ [—L/2, L/2]. Radius 
fluctuations are governed by the stochastic Ginzburg- 
Landau equation where V(cj>) arises from the zero- 
temperature electron-shell potential of Fig. H3 

Details of the calculations appear in l35| ; here we sim- 
ply summarize the results, which appear in Table |U 

The lifetimes tabulated for sodium nanowires in Ta- 
ble [I] exhibit a rapid decrease in the temperature interval 
between 75K and 100K. These lifetimes correlate well 
with the observed temperature dependence of conduc- 
tance histograms for sodium nanowires [32l IsHL . A 
comparison of the lifetimes of sodium and gold nanowires 
listed in Table [I] indicates that gold nanowires are much 
more stable. In our model this arises from the differ- 
ence in surface tension: aj^ u = 5.9 ciNa, and is consis- 
tent with the observation that noble metal nanowires are 
much more stable than alkali metal nanowires. 

There is also an important prediction contained in Ta- 
ble HJ namely that nanowire lifetimes, which exhibit sig- 



FIG. 6: Electron-shell potential V^heii(-R) at three tempera- 
tures, computed from the free-electron model of Ref. |3q|. 

nificant variations from one conductance plateau to an- 
other, do not vary systematically as a function of radius. 
It can be seen from Table [I] that the activation barriers 
vary by only about 30% from one plateau to another, and 
that a wire with a conductance of 96Go has essentially 
the same lifetime as that with a conductance of 3Gq. In 
this sense, the activation barrier variation exhibits uni- 
versal mesoscopic fluctuations: in any conductance inter- 
val, there are very short-lived wires (not shown in Table 
P) with very small activation barriers, while the longest- 
lived wires have activation barriers of a universal size: 

/t2 \ 1/2 

< \E OQ < 0.7 , (25) 

\m e J 

The derivation of l|25|l will be presented elsewhere |37j . 

At present, lifetimes of wires shorter than L c have not 
been systematically studied. We expect that over the 
next several years the technology will improve to the 
point where this will become possible, and our predic- 
tions for a transition in both activation barriers (from 
barriers essentially independent of length to those vary- 
ing linearly with length, as in Fig. and prefactors 
(Fig. can be tested. 

B. Magnetic reversal in nanomagnets 

The dynamics of magnetization reversal in submicron- 
sized, single-domain particles and thin films is important 
for information storage and other magnetoelectronic ap- 
plications. This problem can be treated with the methods 
used throughout this paper, but with a more complicated 
equation of motion than Q. The magnetization dy- 
namics is governed by the Landau-Lifschitz-Gilbert equa- 
tion |3Sj perturbed by weak spatiotemporal noise: 

d t M = - 7 [M x H cff ] + (a/M )[M x d t M] , (26) 
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TABLE I: The lifetime r (in seconds) for various cylindrical sodium and gold nanowires at temperatures from 75K to 200K. 
Here G is the electrical conductance of the wire in units of Go = 2e 2 /h, L c is the critical length, and AE X is the activation 
energy for an infinitely long wire. From Ref. pjEf . 



G 

[Go] 


Na 


Au 


M 


[meV] 


r[B] 


L c 


[meV] 


r [s] 


75 K 


100 K 


150 K 


[nm] 


100 K 


150 K 


200 K 


3 


1.4 


210 


30 


9 x lO" 3 


3 x 10 -6 


1.5 


470 


10 11 


2 x 10 3 


0.2 


6 


2.6 


170 


0.06 


9 x 10~ 5 


io- 7 


3.0 


310 


10 3 


7 x 10 -3 


2 x 10" 5 


17 


3.1 


230 


500 


0.08 


10 -5 


3.4 


470 


9 x 10 10 


10 3 


0.2 


23 


3.7 


190 


4 


2 x 10" 3 


10~ 6 


4.1 


390 


10 7 


3 


2 x 10" 3 


42 


4.3 


210 


50 


0.01 


4 x 10~~ 6 


4.8 


440 


3 x 10 9 


100 


0.03 


51 


4.5 


150 


7 x 10" 3 


2 x 10" 5 


5 x 10" 8 


4.9 


320 


2 x 10 3 


0.01 


3 x 10" 5 


96 


5.8 


200 


5 


2 x 10~ 3 


10" 6 


6.3 


440 


8 x 10 9 


300 


0.05 



where Mq is the magnitude of the magnetization M, a 
the damping constant, and 7 > the gyromagnetic ratio. 
The effective field H c ff = —SE/SM. is the variational 
derivative of the total energy E, which is given by (with 
free space permeability /io = 1): 

E[M(x)] = A 2 f d 3 ir|VM| 2 + f d 3 x\\7U\ 2 

- f d 3 xH c -M, (27) 
Jn 

where f2 is the region occupied by the ferromagnet, A is 
the exchange length, and U (defined over all space) satis- 
fies V ■ (WU + M) = 0. The first term on the RHS of J23 
is the bending energy, the second the magnetostatic en- 
ergy, and the last the Zeeman energy. We take H c || 6, 
using cylindrical coordinates (f , 9, z). The magnetostatic 
energy is nonlocal and gives rise to shape anisotropics 
(for 'soft' magnetic materials, like fee Co or permalloy, 
crystalline anisotropies are negligible). The out-of-plane 
anisotropy energy is especially strong, and forces M to 
lie in the plane. 

The presence of the nonlocal magnetostatic term com- 
plicates analysis. However, the quasi-2_D nature of the 
problem allows a significant simplification, as shown by 
Kohn and Slastikov ,39] . Their asymptotic scaling anal- 
ysis applies to the present problem when both the aspect 
ratio t/R (ring thickness divided by ring mean radius), 
and A/i?, are sufficiently small. Because of the high en- 
ergy cost of variations in Mq, and because the geometry 
under consideration admits nonsingular solutions for the 
vector field M, Mq can be taken to be fixed. The result 
for the total magnetic energy is then [iof : 



ds 



(^) 2 +: 
OS 



2h cos ( 



(28) 



where £ is related to ring circumference, h is the exter- 
nal magnetic field magnitude, and <f> is the local angle 
between the in-plane magnetization vector and the local 



field direction. In (|28|l energy, length, and field are all 
dimensionless, normalized by a corresponding character- 
istic quantity determined by various ring parameters. 

Eq. I|26l) and the variational equation H e ff = —5E/51sA 
yield a nonlinear differential equation that must be sat- 
isfied by any time-independent solution: 

d 2 (j)/ds 2 = sin cos + ft, sin . (29) 

There are three 'constant' ((f) independent of 9) but 
nonuniform (m = M/Mo varies with position) solutions 
for < h < 1: the stable state <j> = (m = 9); the 
metastable state <fi = tt (m = — 9), and a pair of degener- 
ate unstable states <fi = cos _1 (— h), which constitute the 
saddle for a range of (£,h). The 4> — 0, tt solutions are 
degenerate when h = 0, and the <p = n solution becomes 
unstable at h = 1. 

We have also found a nonconstant 'bounce' solution 
of (|29|l . which is the saddle for the remaining range of 
(£,h). It is 

sn(lZ\m) ' 



(s, m) — 2 cot 



dn 



s - s 



cn(1Z\m) 



(30) 



where dn(-|m), sn(-|m), and cn(-jm) are the Jacobi ellip- 
tic functions with parameter m |28| . Sq is a constant, and 
1Z and S are given by 

sn 2 (7£|m) = l/m-h/2 



(I /2m) ^fm 2 h 2 + 4(1 - m) 

™2 



S 2 = 



2 - (m + ^Jm 2 h 2 + 4(1 - m)) 



(31) 
(32) 



Imposition of the periodic boundary condition yields a 
relation between £ and m: 



2K(m)d. 



(33) 



As m — * 0, dn(a;|0) — > 1, and the bounce solution reduces 
to the constant state = cos _1 (— h). The critical length 
and field where this occurs are related by 



ttS c = 2^/Vl - h 2 c . 



(34) 
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FIG. 7: Total switching rate (in units of s _1 ) vs. inverse 
temperature j3 (in units of °K~ X ), at fields of (a) 52.5 mT 
(nonconstant saddle) and (b) 72.5 mT (constant saddle). Pa- 
rameters used are k = .Of, / = .05, R = 200 nm, Ri = 180 
nm, R 2 = 220 nm, Mo = 8 x 10 5 A/m (permalloy), a = .01, 
and 7 = 1.7 x 10 T~ 1 s . Deviation of low-field switching 
rate in (a) from dashed line signals non-Arrhenius behavior. 
(From Ref. H@|.) 



What is interesting here is that the transition is gov- 
erned by two parameters: not only the length, but also 
the magnitude of the externally applied magnetic field. 
While the former cannot be varied continuously, the lat- 
ter can, allowing for the first time a detailed experimental 
probe of the transition. Fig. |7| shows theoretical predic- 
tions of the magnetization switching rate at two different 
field strengths for a ring of fixed circumference. 



VI. CROSSOVER FROM THERMAL 
ACTIVATION TO QUANTUM TUNNELING 

Only purely classical activation processes have been 
considered so far, but the transition in activation behav- 
ior as interval length is increased has an interesting par- 
allel with the crossover from classical activation to quan- 
tum tunneling as temperature is lowered. The correspon- 
dence can be made explicit by mapping system length L 
in the former case to temperature T in the latter. 

The partition function Z at inverse temperature /3 can 
be written as a Euclidean path integral over trajectories 
q(r), where r = it is imaginary time, with each path 
weighted by its Euclidean action Se |4l| : 



Z = J V[q(r)}exp{-S E [q(r)]/n}. 



(35) 



The integral in (|35[1 runs over all paths periodic in imag- 
inary time, with period h/3. 

Eq. I|35|) is derived from the usual definition of the par- 
tition function Z = Tr {exp[— [3H.]}, where H. includes 
both the system and its environment. Consequently, 
even if one is dealing with the tunneling of only a sin- 
gle degree of freedom q(r), the effects of friction due to 
its coupling with the environment must be included. The 
proper treatment of the effects of damping on quantum 
mechanical tunneling have been considered by a number 



of authors; see, for example, 0,0,0,0- Although 
friction strongly affects the tunneling rate quantitatively, 
it will not be included here; this is because our only aim is 
to present the connection between the transition in ther- 
mal activation of (infinite-dimensional) classical fields as 
length is varied, and the crossover from classical activa- 
tion to quantum tunneling in (one-dimensional) systems 
as temperature is varied. 

An early treatment (that ignored dissipation) was 
given by Goldanskii [12\. Setting the classical, 
temperature-dependent Arrhenius factor AE/ksT equal 
to the zero-temperature quantum tunneling rate through 
a parabolic barrier, he noted that the characteristic tem- 
perature To for the quantum tunneling <-> classical acti- 
vation crossover was 



T = hLO c /(2TTk B ) ■ 



(36) 



where u c is the characteristic frequency of the locally 
quadratic barrier. Although Goldanskii's approach led 
only to an estimate, his formula H3bfl was quite accu- 
rate (in the absence of dissi pati on) , as the following 
more detailed analysis Qa, Ir7l llS| will demonstrate (see 

also 000 0,0 tab- 

In what follows, we will use the asymmetric potential 
given by JSJ. This allows us to consider only incoherent 
tunneling processes: once the particle tunnels through 
the barrier, it escapes for good. A double-well potential 
such as (0} requires consideration of dissipative quantum 
coherence effects, for which a real-time functional integral 
approach is better suited (see, for example, |42ll43| ). The 
Euclidean action for a particle of mass M is therefore 



Se[q(t)] 



pn/2 



dr 



ph/2 

2 



{hlq 2 (r) (37) 

~i^q{ T ) 2 ~ ~^q{ T ) 3 | + (frictional terms) 

where q denotes a derivative with respect to imaginary 
time r. In the low-friction limit, extremal paths satisfy 



M 



q(r) + ulq(j) - \q{rf 



0. 



(38) 



with periodic boundary conditions q{— flh/2) = q((3H/2). 

Three extremal solutions are physically relevant. The 
first is the uniform q(r) — solution, which corresponds 
simply to the stable state at all temperatures. Of the 
other two, one is uniform: 



?c,+(t) = ujI/\ 
and the other is a nonuniform bounce: 

q c ,-(T,m) = -^-v(m) 2 dn 2 ^ ±-L( T - r ) 



(39) 



2A 



1 



m)v(m) 



(40) 



where as before < m < 1 and v(m) = (1 — j7i + to 2 )~ 1//4 . 
Imposition of the periodic boundary condition relates the 
parameter m to the inverse temperature: 



/? = 4K.(m)/hujov(rn) 



(41) 
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the high-temperature limit and the quantum tunneling 
formula is recovered in the zero-temperature limit : 



0.5 1 
(a) 



FIG. 8: (a) The Euclidean action SE[g c (r)] at all temper- 
atures, (b) The logarithm of the leading-order escape rate, 
shown in an Arrhenius-style plot. In both graphs, M = loo = 
A = H = ka = 1, and the dot indicates the crossover. 



The constant tq in l|4(J|l indicates the translational sym- 
metry in imaginary time arising from the periodic bound- 
ary conditions, and corresponds to a zero mode as de- 
scribed in Sec. IIIII 

We can now make explicit the mapping to the thermal 
activation of classical fields. Here the inverse tempera- 
ture f3 plays the same role as system length L in Sec. II VI 
In particular, high temperature corresponds to the small 
length regime. For /3 smaller than some (3 C , we would 
therefore expect the constant solution q c + to be the sad- 
dle. And indeed it is. The corresponding action is 

Ssfec+e-)] = Mtu%l3h/(6\ 2 ) = (3HAV , (42) 

where AV — AIluq/ (6X 2 ) is the potential energy differ- 
ence between the potential barrier top and bottom. Con- 
sequently, exp{— SE[q c ,+ (T)]/h} = exp[—/3AV], and the 
classical Arrhenius factor is recovered. 

The bounce solution Se[iIc,-(t)] is the saddle at large 
(3. As in Sec. IIVI the transition temperature T c is given 
by gTJ when m -> 0+. We find that T c = fiw /(27rfc B ), 
exactly that found by Goldanskii (Eq. i|3t)|) h using much 
simpler arguments. (For the potential used in (|37[) . the 
curvatures at the top and bottom of the well are equal.) 

The action in the low-temperature region 

S E [q c ^} = 3Afco 5 /(5A 2 (l - m + m 2 ) 1 ^) 
x [2E(m) - [(2 - m)(l - m)/(l - m + m 2 )]K(m 



+{Mj3hu)l/ '12A 2 ) 1 - 3(2 - m)/2\/l - m + m 2 

+ (2 - m) 3 /2(l - m + m 2 ) 3/2 j (43) 

is, as in the classical field case (cf. Fig. |2J|, continuous 
and differentiable at all temperatures; but also as before, 
its second derivative is discontinuous. The action, along 
with the leading-order escape rate, is shown in Fig. |S] 

The well-known zero-temperature tunneling rate is re- 
covered from (|43|l in the limit m — > 1. Summarizing, we 
find that the classical Arrhenius formula is recovered in 



exp 



-S E {q c ,-)/ft 



exp[-AV/k B T] T -> oo 



(44) 



exp[-36AI//5^ ] T -> 
while in a narrow region about Tq both contribute. 

VII. DISCUSSION 

We have demonstrated that in problems involving 
noise-induced escape of a classical field over a barrier, a 
type of second-order phase transition, with what appears 
to be attendant critical phenomena (such as power-law 
divergence of the rate prefactor) occurs as one or more 
external parameters (length of the interval on which the 
field is defined, external magnetic field if relevant, and 
so on) is varied. We have also shown that there is a 
mathematical mapping of this transition to the classical 
activation <-» quantum tunneling transition for a particle 
escaping a simple potential well. The mapping here in- 
volves identifying interval length (and/or magnetic field, 
if appropriate) in the classical field case to temperature 
in the classical <-> quantum case. To avoid confusion, it 
should be remembered that the asymptotically small pa- 
rameter is noise strength (typically, but not necessarily, 
temperature) in the former problem, and Planck's con- 
stant — not temperature — in the latter. 

In this section we will address two questions that im- 
mediately come to mind: To what extent can the transi- 
tions discussed be considered 'real' second-order phase 
transitions exhibiting critical phenomena; e.g., in the 
sense that the fluctuations driving the transition occur 
on all scales? Secondly, how deep is the correspon- 
dence between the classical field transition and the quan- 
tum <-* classical crossover O E3, EE E3, EE El 

m? 



A. Is the phase transition 'real' ? 

The short answer is: in a mathematically asymptotic 
sense yes, but from a strictly (and more physically rele- 
vant) critical-phenomena-oriented viewpoint, no. More- 
over, the observation of rate prefactor divergence depends 
crucially on the order in which relevant parameters (tem- 
perature, system length, and so on) are varied. There are 
clearly different saddle configurations, and therefore qual- 
itatively different activation behaviors, on either side of 
the transition. But the question we are focusing on here 
is: what is happening very close to the critical length- 
scale? This was extensively discussed in and that 
discussion will be expanded here. 

Naively, the transition appears second-order in several 
respects: the saddle solution is continuous (and even bi- 
furcates in symmetric models such as (@J)), and the action 
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is continuous and differentiable at the transition point 
but has a discontinuous second derivative there. But 
perhaps most compelling, from the point of view of criti- 
cal phenomena, is the apparent power-law 'divergence' of 
the rate prefactor shown in Fig. [5] We therefore examine 
this in more detail, first asking what it even means for 
the prefactor to 'diverge'. 

Of course, at no lengthscale is the actual prefactor 
infinite. Consider the analysis of the noisy symmetric 
Ginzburg-Landau model given in Sec. II VI It is important 
to recall that the analysis of the escape rate is, strictly 
speaking, valid only in an asymptotic sense as T — > 0: 
our results apply only to temperatures sufficiently low so 
that the escape rate is small. 

What the formal divergence of the prefactor does mean 
is that the escape behavior becomes increasingly anoma- 
lous as L c is approached, and that it is non-Arrhenius 
exactly at L c , where for all T — > the prefactor is 
temperature-dependent, scaling as a negative power of T. 
In the region close to L c , the rate prefactor Tq is anoma- 
lously large, but still finite. The formula 1|12|) — from 
which the prefactor shown in Fig. [S] was computed — is 
valid only for T sufficiently small so that the contribu- 
tions from the quadratic fluctuations about the relevant 
extremal state of H[<f>] dominate the action. So as long 
as all eigenvalues of A t are nonzero (excluding, if transla- 
tional symmetry is present, the zero mode which may be 
extracted), Eq. II 21 ) applies, but only within a tempera- 
ture region driven to zero as L — > L c . The diminishing 
size of this region as L — > L c is controlled by the rate of 
vanishing of the eigenvalue(s) of smallest magnitude. 

The implication is that the prefactor behavior depends 
on whether T — > at fixed L near L c , or whether L in- 
creases (say) through L c at fixed low T. In the former 
case, one will recover Fig. [5] If instead one fixes tempera- 
ture at some small but nonzero value, one should observe 
first a rising prefactor as L approaches L c , but at some 
L (depending on T through a type of 'Ginzburg crite- 
rion'; cf. Fig. 9 in 11]), the prefactor crosses over to a 
non-Arrhenius (temperature-dependent) form. As L con- 
tinues to increase on the other side of L c , the sequence 
of events is reversed. 

The procedure of fixing T and varying L is in many 
instances the more physical one. In this case a correct 
analysis needs to include higher-order (than quadratic) 
fluctuations about the transition state, as was done, e.g., 
in [l7|. The next higher-order terms will (unless pre- 
vented by symmetry) be nonzero. The behavior remains 
anomalous, however. Suppose one starts to vary L at 
fixed T, but as soon as non-Arrhenius behavior is en- 
countered, one fixes L and then starts to lower T. A 
subsequent transition from non-Arrhenius to Arrhcnius 
behavior will again be encountered: as T is lowered, the 
prefactor Tq rises until it reaches the value shown in 
Fig. [5J it remains constant thereafter. 

Lying behind this description is the relative magni- 
tudes of the thermal energy and that due to quadratic 
fluctuations in the direction of the eigenmode r\\ with 



vanishing eigenvalue \\. Non-Arrhenius behavior should 
be seen at 'intermediate' temperatures, where thermal 
energy is large compared to that due to quadratic fluc- 
tuations along the r\\ direction, but small compared to 
that arising from higher-order fluctuations. 'Intermedi- 
ate' here depends on L, whose closeness to L c determines 
the magnitude of Ai. Arrhenius behavior reappears in 
the 'low' temperature region where the thermal energy is 
small compared even to that due to the quadratic fluc- 
tuations along 771 . And because Ai — > as L — > L c , the 
crossover from non-Arrhenius to Arrhcnius behavior oc- 
curs at an increasingly lower temperature. Exactly at 
L c , where Ai is exactly zero, the escape behavior is non- 
Arrhenius for all T — > 0. 

There's an alternative way of describing the situa- 
tion: when viewed on 'normal' fluctuation lengthscales of 
0(T 1 / 2 ), field fluctuations along the eigenmode direction 
771 appear to be diverging as L — > L Cl accompanied by 
anomalous transition behavior. However, when viewed 
on an 'anomalous' lengthscale ((^(T 1 / 3 ) or (^(T 1 / 4 ) or 
depending on the form of the potential), those fluctu- 
ations remain finite, and one would presumably observe 
a rounded maximum of the prefactor (when scaled by the 
appropriate temperature-dependent factor) at L c . 

The situation is therefore fairly subtle; care must be 
used to describe the 'transition' in appropriate terms. 
Moreover, in perhaps the most important respect, the 
transition fails a central test of criticality — that of 
the disappearance of a characteristic fluctuational length- 
scale. In our problem, even at L c there is a well-defined 
characteristic lengthscale, albeit an anomalous one. 

This is strikingly different from a similar-loo king 
(on the surface) transition described elsewhere 0, l45l 
|46| . This transition in the activation behavior of non- 
equilibrium systems (e.g., where detailed balance is ab- 
sent in the zero-noise dynamics) occurs as the result of 
singularities developing m the action. These singu- 
larities lead in turn to the appearance of caustics in the 
pattern of fluctuational trajectories, as in geometric op- 
tics. Here the transition occurs as a parameter in the 
zero-noise dynamics is varied. 

In these systems, fluctuations about the optimal es- 
cape trajectory do occur at all scales at the critical 
point, so one cannot simply 'cure' the divergence by in- 
cluding higher-order terms as before. A de novo scaling 
theory |46t I47I ] is required. This theory results in an 
array of nontrivial critical exponents that obeying scal- 
ing relations and characterize the divergence (or vanish- 
ing) of relevant physical quantities (including, but not 
limited to, the rate prefactor). Consequently, these non- 
equilibrium systems can justifiably be said to exhibit true 
critical behavior at the transition point. 

The discussion so far has focused entirely on second- 
order transitions; what about first-order? This possibil- 
ity has been discussed by several authors 0, |2(], 0] • 
In slightly more complicated classical field theories per- 
turbed by spatiotemporal noise, such as a sixth-degree 
Ginzburg-Landau model, the nonconstant saddle branch 
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of the energy functional TL[-\ can in principle cross the 
constant saddle branch at a nonzero angle. This should 
give rise to a first-order transition. So, in the phase plane 
of these models, the second-order transition point (in the 
limited sense discussed in this section) is presumably the 
endpoint of a first-order transition curve. 

B. Is the quantum tunneling «-» classical activation 
crossover for a single degree of freedom identical to 
the transition in activation behavior of a classical 
field on finite intervals? 

Mathematically, yes, under the following mapping: 

Classical <-> Quantum Classical <j> 4 

Small parameter h T 

Tunable parameter T L 

Periodic in (3h L 

Here we chose a stochastic Ginzburg-Landau 4 model 
for specificity, but one can substitute any of the other 
classical field theories discussed. For the noisy magneti- 
zation dynamics discussed in Sec. IV Bl it should be re- 
membered that the transition can occur as either length 
or field is varied. One also need not choose periodic 
boundary conditions for classical Ginzburg-Landau field 
theories; the transition will occur in similar fashion (but 
with fairly minor differences as described in Sec. UVfl with 
other types of boundary conditions. On the other hand, 
one is constrained to use boundary conditions periodic in 
(3h in the quantum «-> classical problem. 

This mapping is realized when the quantum <-> classical 
transition problem is set up using a Euclidean imaginary- 
time functional integral formulation, in a potential where 
incoherent tunneling dominates. The classical problem is 
approached similarly as a real-space path integral in the 
limit of weak spatiotemporal noise. 

But there is a significant physical difference (leaving 
aside the obvious ones) between the two classes of prob- 
lems. In the previous section, we were able to consider 
— for the classical problem — either holding L fixed and 
lowering T, or holding T fixed and moving L through 
L c . In the first case one can, in principle, recover the 
prefactor divergence shown in Fig. \5\ while in the sec- 
ond one should observe a sequence of Arrhenius — ► non- 
Arrhenius — > Arrhenius transition behaviors, with the 



width of the non- Arrhenius region vanishing as T — > 0. 
'Non- Arrhenius' means, as usual, a prefactor T not de- 
termined by l|12p. and thereby usually aquiring a tem- 
perature dependence. In order to see this dependence, 
one would need to repeat the procedure (fix T, move L 
through the transition) at several different temperatures. 

But in the quantum case, there is no such freedom: 
here the small parameter, h, cannot be varied. So a pref- 
actor divergence such as Fig. JSJ) cannot be observed even 
in principle. 

In a similar vein, the zero mode arising from the use of 
periodic boundary conditions will differ in the two cases. 
This mode arises on the nonconstant ('bounce') side of 
the transition due to translational symmetry — the in- 
stanton (or domain wall) separating the two stable solu- 
tions can arise anywhere in space (or imaginary time, in 
the quantum *-y classical case). This leads in turn to a 
T-dependent (respectively, ^.-dependent) prefactor for all 
L > L c (respectively, T <T C ). 

It should be kept in mind that we have made the im- 
plicit assumption in the classical field case that temper- 
ature can be lowered to arbitrarily small values. For 
many systems this may not be the case: either tran- 
sition rates will become immeasurably low, or the sys- 
tem itself might undergo a (real) phase change, or new 
physics may enter in some other way. Nevertheless, the 
ability to vary temperature does distinguish the classical 
field transition from the quantum <-> classical one. In 
particular, the magnetic reversal problem (cf. Sec. IV B|) 
presents us with the exciting possibility of continuously 
tuning through the transition by varying an external 
magnetic field. With the possibility this presents of ex- 
perimentally studying activation closer to the transition 
than might otherwise be achieved, we might in this way 
not only advance our understanding of stochastic reversal 
in nanoscale magnets, but also uncover new information 
about the quantum <-* classical transition that the con- 
stancy of h would otherwise prevent us from obtaining. 
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